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We take advantage of the wealth of rotation measures data contained in the NRAO VLA Sky 
Survey catalogue to derive new, statistically robust, upper limits on the strength of extragalactic 
magnetic fields. We simulate the extragalactic magnetic field contribution to the rotation measures 
for a given field strength and correlation length, by assuming that the electron density follows the 
distribution of Lyman-a? clouds. Based on the observation that rotation measures from distant radio 
sources do not exhibit any trend with redshift, while the extragalactic contribution instead grows 
with distance, we constrain fields with Jeans’ length coherence length to be below 1.7 nG at the 20 - 
level, and fields coherent across the entire observable Universe below 0.65 nG. These limits do not 
depend on the particular origin of these cosmological fields. 


MOTIVATION 

Is the Universe permeated by an all-encompassing 
magnetic field (MF)? MFs are already observed in nearly 
all types of structures, from planets to galaxies to clusters 
of galaxies CHI], but cosmological magnetic fields still 
remain elusive. Nonetheless, understanding their char¬ 
acteristics has cardinal relevance in many fields in as- 
troparticle physics and cosmology: propagation of ultra- 
high energy cosmic rays (UHECRs), structure formation, 
early and very early Universe models, physics beyond the 
Standard Model, radio-astronomy, and so on, see m- 

One way cosmological MFs manifest themselves is 
by rotating the plane of polarisation of electromagnetic 
waves propagating from far away sources to the Earth. 
The main idea of this work stems from a very simple 
observation we made in our previous paper [7]: Fara¬ 
day rotation measures (RMs) of distant objects do not 
show any evolution with redshift. However, if there is an 
all-pervading, extragalactic MF (egMF) we do expect a 
quite pronounced change in the distribution of RMs with 
distance. We can thus limit the strength of such a field 
by comparing simulated and observed RMs distributions. 

Currently, the strongest upper limits on the strength, 
B , of present-day egMFs come from microwave back¬ 
ground observations [5] (see also |5J [ID]) and read B < 
2.8 nG for a coherence length l c — 1 Mpc. In the special 
case of a scale-invariant MF spectrum, these limits could 
be further lowered to 0.9 nG, if one accounts for their 
impact on the ionisation history of the Universe. Notice 
that these limits apply only to primordial MFs, i.e., fields 
generated in the very early Universe, while, for instance, 
cosmological MFs could be generated at later stages by 
various astrophysical mechanisms (e.g. m). 

The limits coming from the analysis of RM data are 
less restrictive: according to [14.] (see also [15]), fields 
correlated on Mpc scales are bound to have B < 6 nG. A 
comprehensive overview of present and earlier constraints 


coming from RM observations can be found in [3] . 

We devise a method for extracting the unknown er¬ 
rors from the data itself, and thanks to this new method 
and the new available RM data, we can improve on these 
limits by five times, and include a full treatment of their 
statistical significance: we find that egMFs with coher¬ 
ence lengths of about 1 Mpc and strengths above 1.7 nG 
(1.2 nG using only low-luminosity (Ip) sources) are in¬ 
compatible with current RM observations at the 2er level; 
this limit becomes 0.65 nG (0.50 nG) if the egMF is co¬ 
herent across the entire Universe. 

METHOD 

Observations. The plane of polarisation of a linearly 
polarised electromagnetic wave which moves through a 
magnetised plasma rotates by an angle (p proportional to 
the square of the wavelength X: ip = RMA 2 , where 


n e (z)B\\(z) 

d l(z) 

Id (1 + ^) 2 

dz 


Here n e is the density of free electrons measured in cm' 3 , 
U|| is the component of the MF (in pG) parallel to the 
line of sight l(z), and D is the distance to the source in 
kpc; here and everywhere RM is measured in rad/m 2 . 
Notice that in order to gain some knowledge about the 
MF, some independent estimates of n e are required. 

The largest set of RM of extragalactic sources to date 
was compiled in [T6] from the NRAO VLA Sky Survey 
(NVSS) data E]. The total number of observed sources 
was 37,543, of which 4002 have known redshifts m- 
From this set, we accepted only sources with galactic 
latitude \b\ > 20°. 

Generically, we can split the observed RM as RM 0 b s = 
RMg a i + RRM, where the first term is the contribution of 
the regular MF of the Milky Way, and the second term 
stands for “residual RM”, and encodes all other sources 
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Zb 

0.15 

0.35 

0.7 

1.3 

1.65 

1.95 

2.25 

2.6 

5 

N s 

418 

418 

501 

677 

291 

137 

76 

50 

25 

(|RRM|) 

16.2 

15.3 

15.9 

16.6 

15.4 

15.8 

16.2 

13.9 

16.3 


TABLE I. Upper bin redshift boundaries Zb, numbers of 
sources in the bin N s , and their averages (|RRM|), for the 
low-power set. 


of RM once the local (regular) MF is subtracted: RM 
instrinsic to the source, measurement errors, turbulent 
galactic MF, and egMF (RRMx), see [T9] for an estimate 
of these components. 

In order to disentangle redshift-dependent effects from 
those pertaining to the sources themselves, we computed 
their luminosity using, where available, the most recent 
spectral indices from [20]; since not all sources have mea¬ 
sured spectral indices, this reduced the total number of 
objects to 3053. This set was split into two using a 
luminosity threshold of Ti. 4 GHz = 10 27 8 W Hz -1 : the 
(Ip) group counts 2593 sources, while the high luminos¬ 
ity one includes 460 of them — notice that this particular 
choice of luminosity cut does not affect our results sig¬ 
nificantly, see [7]: we have also checked this statement 
directly. When the RM ga i contribution is subtracted 
from RM ohs m, one can see that RRMs of sources of 
lower power have an evolution with redshift consistent 
with zero. High power sources are systematically shifted 
towards higher RRM, but, with the data available to¬ 
day, it was not possible to determine whether this set 
evolves with z. Table |T] summarises this observation for 
the choices of luminosity cut and redshift binning which 
were made in [7]. In this work we employ both the full 
and the low luminosity sets; we also tested that our re¬ 
sults are anyhow robust against the selection of the cut. 


Simulations. When the extragalactic medium is per¬ 
meated by a MF, this leaves an imprint on the observed 
RMs, see for example m the RRMx systemati¬ 
cally grow with redshift due to their accumulation along 
the line of sight. The featureless behaviour which we ob¬ 
served instead is thus incompatible with what is expected 
if an egMF were present. By simulating the effects of this 
hypothesised egMF we can thence constrain its strength. 

In order to build a model for RRMx we need to specify 
the properties of the egMF and the electron density n e . 
Following HUES], the latter is assumed to be well de¬ 
scribed by the observed Lyman-a (Lya) forest distribu¬ 
tion of neutral hydrogen absorption lines. In particular, 
we take the analytical approximation which was given 
in HU EH, which is a standard log-normal distribution 
for the electron overdensity , S e , with scale parameter 


a e (z ) = 0.08 + 
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4.21 | 1.44 

(1 + 2) 2 + (1 + 2) 3 ’ 


and location p e (z) = —ct 2 (z)/2: 
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This distribution is accurate for fluctuations at the 
Jeans length scale Aj(z) ~ 2.3(1 + z)~ 3 / 2 Mpc (23] [251 - 
[27] (we adopt Hq = 71 km/s/ Mpc as the Hubble 
parameter today and f 1m = 0.27 as the total matter 
density fraction) [25]. The actual electron density is fi¬ 
nally expressed as n e (z) = n e (0)(l + <5 e )(1 + z) 3 , with 
n e ( 0) = 1.8 x 1(U 7 cm -3 . 

The MF is characterised by the strength B which we 
wish to constrain, and the power spectrum whose shape is 
in principle not known. However, for most cases it can be 
checked that the effect of the MF on the observed RM is 
dominated by a single scale, the coherence length l c which 
we treat as a free parameter; this scale is not necessarily a 
physical scale, but a useful practical tool for the analysis: 
our results are easily rescaled for/adapted to a specific 
MF model (see [29] )• In our simulations we test values for 
l c from a tenth of the Jeans length A j and the Hubble size 
1 /Hq. Since the conductivity of the Universe is extremely 
large, and since diffusion of the MF is inefficient at scales 
much larger than 1 AU, we can safely assume that the MF 
is frozen into the plasma (neglecting non-linear effects), 
see for instance 00; thus, for spherical overdensities, it 
will scale accordingly as B(z) = B [n e (z)/n e (0)] 2 / 3 . 

Practically, we generate a large number of lines of sight 
in steps of A j up to some given redshift, and we collect 
RRMx from each step. The electron density is generated 
at each step sampling the distribution © ; to simulate 
the randomness of the MF orientation, we recalculate its 
amplitude each time the distance travelled equals a mul¬ 
tiple of the correlation length, drawing from a uniform 
[—1,1] distribution M- 

With these prescriptions, we have first obtained the 
expected theoretical egMF-induced |RRMx| evolution 
curves with redshift, by averaging 10 5 lines of sight out 
to redshift 2 = 5: the result is given in Fig. □for a bench¬ 
mark MF reference value of B = 1 nG: the rapid increase 
of |RRMx| with redshift is evident. Moreover, there is a 
clear transition, due to the overall redshift-dependence 
of | RRMx|, from lower to higher redshifts, roughly lo¬ 
calised between z = 0.5 and 2 = 1: for l c = 1/Hq the 
growth of RM is significantly damped at high redshifts, 
whereas if l c = A j the curve becomes essentially flat. 
This transition is also observed in the different shapes of 
the |RRMx| distributions at different redshifts, and is in 
agreement with the earlier results of HUES- 

In our simulations we can generate only the contribu¬ 
tion from the egMF. However, as we mentioned before, 
when comparing with the data one needs to account also 
for the turbulent Galactic fields, measurement errors, and 
intrinsic RMs. The two first contributions are dominant 
and do not depend on redshift; the third one does but it is 










3 



z 


FIG. 1. Theoretical egMF-induced |RRMx| evolution with 
redshift for l c = A j (red, solid), and l c = 1/Ho (blue, dashed), 
averaged for 10 5 lines of sight. 

subdominant — it is precisely this feature, which depend 
on redshift in a very specific way, that allows us to make 
this statement, since it is not seen in the data. We can 
then exploit the data at low redshifts to obtain infor¬ 
mation about these redshift-independent contributions, 
also since at low redshifts they are going to be more im¬ 
portant, and then utilise the high-redshift portion of the 
data to compare with the simulated distributions for the 
same set of sources, with the egMF included. For another 
approach to extract the extragalactic RM piece see m- 

We thus split the data in three redshift bands: for sim¬ 
plicity, and because it well matches what we would infer 
from Fig. |T| and more consistently by directly compar¬ 
ing the underlying distributions, we take the low redshift 
band to correspond to the bins 1 and 2 of Table [1] that is, 
2 = [0,0.35], the high redshift band from z > 1.3 (that 
is, bins 5 to 9), and a transition band — which we do not 
use — corresponding to bins 3 and 4, or z = [0.35,1.3]. 
The three sets contain 836 (836 in the (Ip) set), 1254 
(1178), and 936 (579) sources, respectively. Reducing the 
size of the intermediate band does not strengthen our 
constraints but marginally, despite the fact that many 
sources are found in this redshift range, because the 
increase in number of sources is counterbalanced by a 
less pronounced contribution of egMF compared to other 
sources of RRM. 

We finally build the needed simulated distribution of 
|RRM| at high-z as follows. First, we randomly pick one 
RRM from the low-z set: this serves as our estimation 
of the ^-independent contribution which we can not sim¬ 
ulate in our model (and for which a specific value we 
do not want to adopt a priori)', this is possible because 
in this low-z bin we have that |RRMx| -C |RRM|. In 
order to obtain the total |RRM|, we then generate a sec¬ 
ond batch of RRMx values by simulating 100 lines of 
sight for each of the sources of the high-z set; that is, 


the RRMx of all high-z sources is simulated 100 times 
for each source, out to their actual redshift. This gives 
a total of 96,300 (57,900) lines of sight and correspond¬ 
ing RRMx. These RRMx are generated for our bench¬ 
mark field value of B = 1 nG, and then rescaled for any 
other value of B , simulating the egMF contribution we 
are seeking. Randomly picked values from both batches 
(one each) are then incoherently added (that is, each with 
its own sign) 10 s times to generate the final theoretical 
|RRM| distribution, as a function of the MF strength and 
coherence length, that we can compare with the actual 
data (strictly speaking this procedure is consistent only 
for not too small values of the MF strength in order to 
not introduce ties in the distribution[32]). Notice that 
by definition all errors are included in this procedure as 
we keep the full distribution of |RRM|, not only some of 
its momenta. 


RESULTS 

We compared the two distributions — the data, and 
the theoretical predicted sample — by means of a 
Kolmogorov-Smirnov (KS) test. This test allows us to 
exploit all of the information contained in the distribu¬ 
tions, and is more fit to perform this analysis than a 
simple comparison of means in a given redshift bin, since 
these tend to fluctuate widely due to the underlying log¬ 
normal distribution for n e , which is indeed the major 
source of fluctuations. 

In Fig. [2] we show an example of the PDFs and CDFs 
of the two distributions we are comparing: the data 
((Ip) set) and a simulated |RRM| with B = 3 nG and 
l c = A j, where it is clear that such a field value is strongly 
disfavoured — the two distributions are statistically in¬ 
compatible. 

Fig. I contains the p-values of the KS tests ((Ip) set) 
as a function of B for l c = A j (red, solid) and l c = 1 /Hq 
(blue, dashed) cases. For an egMF with l c = A j we can 
read off the values of B corresponding to p-values of 2a 
and 3 <t as B = 1.2 nG, and B = 1.7 nG, respectively — 
these limits read B = 1.7 nG, and B = 2.2 nG if we keep 
the full data set. For the Universe-wide case the limits 
are somewhat stronger: 2a is already attained at B = 
0.50 nG (B = 0.65 nG for the full set). One may worry 
that in this case however the expected |RRMx| evolution 
with redshift is not constant, see Fig. [T] We have checked 
that the result remains the same if we further split the 
high redshift band into two because neighbouring high 
redshifts generate similar |RRMx| distributions. 

In Fig. [4] we extend the analysis to variable coherence 
lengths; shown here are confidence interval contours ob¬ 
tained from the p-values of the KS tests as a function of 
B (y-axis) and l c = A j (x-axis). We see here how the 
limits can be automatically rescaled for l c < A j, since 
n e fluctuations are dominated by the A j wavelength: 
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FIG. 2. PDF and CDF showing |RRM| for the data (blue, 
dashed) and a simulated |RRM| with B = 3 nG and l c = Xj 
(red, solid). 

the bound on the field strength becomes approximately 
(Xj/lc) 1 ' 2 > 1 times weaker (that is, the allowed egMF 
is (A j/lc) 1 ^ 2 > 1 times stronger). 


FIG. 3. KS tests p-values as a function of B for l c = X j (red, 
solid), and l c = 1/Hq (blue, dashed). 



Coherence length [A j\ 


CONCLUSION 

Extragalactic MF with coherence lengths of about 
1 Mpc can not be stronger than 1.7 nG at the 2cr level 
(1.2 nG for the low-luminosity data), whereas a Universe¬ 
wide egMF is bound to be weaker than about 0.65 nG 
(0.50 nG). These limits are obtained using RM data from 
extragalactic sources, are valid independently of the ori¬ 
gin of these egMF, and for a very large class of egMF 
models. Moreover, these limits are a fivefold improve¬ 
ment over those previously available in the literature, and 
are now more than competitive with microwave back¬ 
ground ones (which however apply only to primordial 
fields). 

This improvement stems from the observation that 
RMs of distant objects do not evolve with redshift. This, 
combined with a much larger and better set of data, 


FIG. 4. Confidence interval contours (bottom panel) obtained 
from the p-values of the KS tests as a function of B (y-axis) 
and l c = A j (x-axis). 


the RM compilation of the NVSS catalogue, enabled a 
much more robust statistical approach to constraining 
the egMF. 

A straightforward application of these limits is to 
UHECRs propagation: if an egMF with B ~ 1 nG and 
l c ~ Aj existed, it would alter the way UHECRs propa¬ 
gate by deflecting them quite significantly. The median 
deflection for a proton primary of even the highest energy, 
10 20 eV, would be around 9 deg when propagating from a 
distance of 200 Mpc. This is a result of a simulation with 
the same Lya distribution for n e and B rather than the 
simplified homogeneous cell model usually implemented. 

Our approach can be directly ported to analyse the 
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data coming from the next generation of radio telescopes, 
in particular the Square Kilometer Array (SKA) [25]. 
The main limitation in the currently available data is 
the low number of sources, particularly those with pre¬ 
cise redshift and spectral index information. The new 
data will significantly overcome this limitation, and al¬ 
low a much more robust determination, and a marked 
strengthening, of our limits; we discuss the sensitivity of 
our new technique in the forthcoming SKA era in a future 
publication. 
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